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Competition between phase coherence and correlation in a mixture of Bose-Einstein 

condensates 
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Two-species hard-core bosons trapped in a three-dimensional isotropic harmonic potential are 
studied with the path-integral quantum Monte Carlo simulation. The double condensates show two 
distinct structures depending on how the external potentials are set. Contrary to the mean-field 
results, we find that the heavier particles form an outer shell under an identical external potential 
whereas the lighter particles form an outer shell under the equal energy spacing condition. Phase 
separations in both the spatial and energy spaces are observed. We provide physical interpretations 
of these phase separations and suggest future experiment to confirm these findings. 
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The rich physics of multi-species boson mixtures has 
attracted much attention since the first realization of the 
coexistence of the double condensates from two differ- 
ent spin states [F = 2,m = 2) and \F = l,m — —1) 
of 87Rb 0. Some interesting properties, such as the 
stability of condensates, phase separation, and symme- 
try breaking, have been investigated both experimen- 
tally and theoretically [1 SE IF, "7, Efl IH [12 . 
It was found that *^Rb in two hyperfine states can form 
two separate condensates with a spatial overlap T| and 
this overlapping region can be adjusted by tuning the in- 
terspecies interaction. Theoretical studies including both 
the Thomas-Fermi approximation and numerical solu- 
tion of the Gross-Pitaevskii equation fl, 0, Q, 0| have 
also found that two overlapping condensates can be sep- 
arated into a core at the trap center and a surrounding 
shell. This phase separation can be achieved by changing 
the interspecies interaction, the total number of particles, 
and the symmetry of the trapping potential. It is sug- 
gested that the system may become unstable for certain 
interaction strength Q , which may lead to a spontaneous 
symmetry breaking spatially 6, 7, 8]. 

Most of the early theoretical efforts mentioned above 
only give the properties of condensates at zero temper- 
ature. In this Letter, we study a two-species boson 
system in a three-dimensional isotropic harmonic trap 
through the path-integral quantum Monte Carlo simu- 
lation, which allows us to probe such a system at finite 
temperature. Here we concentrate on the conditions un- 
der which the phase separation occurs. 

The path-integral quantum Monte Carlo method has 
been successfully applied to the one-species boson sys- 
tems in traps [l^- It provides the exact solution of a 
many-body system within a controllable variance. The 
total energy, density profile, specific heat, condensation 
fraction, and other quantities can all be evaluated [Tsj. 

The system considered here consists of two types of 
particles, marked as species 1 and species 2, respec- 
tively, with total particle numbers A^i and N2 and par- 
ticle masses mi and m2. For a realistic system, these 
two types of particles can be two different elements, two 
different isotopes of the same element, or two different 



hyperfine states of identical atoms. The interactions be- 
tween two like and unlike atoms are characterized by 
the s-wave scattering lengths an, 022, and ai2. In the 
model studied here, an and 022 are the hard-core diam- 
eters of the particles in species 1 and 2, respectively, and 
ai2 = (an + a22)/2. 

The Hamiltonian for the system of = A^i + A^2 par- 
ticles is given by 
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is the Hamiltonian of the noninteracting system of the 
ith species under the spherically symmetric trapping po- 
tential Ui{r) — fcir^/2 ~ nriiLjj'fr'^ /2. The interaction be- 
tween any two particles Vjk is assumed to be a hard-core 
potential. 

The algorithm applied to the one-species systems 
is generalized to study the two-species system here. We 
treat the two-species system as two subsystems, each of 
which contains one species and follows its own statistics. 
The permutations are only performed among the iden- 
tical particles during a simulation. However, these two 
subsystems do not behave independently because of the 
interspecies interaction. The Monte Carlo steps are in- 
fluenced by the interactions among all the particles as 
well as the permutations of the particles in each species. 

In order to calculate the density profile, we divide the 
space into small shells along the radial direction. We 
choose the interval between two shells such that after 
averaging over a large number of Monte Carlo steps on 
the order of 10®, the finite-size error diminishes in com- 
parison with the statistical error. The spatial density 
Pi{r) can then be obtained by counting the number of 
particles from each species in a given shell. The nor- 
malization relation J pi (r) dr = Ni is used to check the 
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convergence. For convenience, we use the dimensionless 
units with h = mi = loi = 1. 

We start the simulation with two species with an iden- 
tical mass, hard-core radius, and total particle number to 
verify that the condensation does occur below a certain 
temperature. The two condensates show different crit- 
ical temperatures as we change the ratios 011/022 and 
TO2/TO1. The stronger the interaction is, the lower the 
critical temperature. The condensation favors the lighter 
particles because a heavier particle system has a lower 
critical temperature, which agrees well with the available 
experimental results . 

Our main interest here is on the structures of the mix- 
tures. We calculate the densities of two condensates at 
various parameters at temperature T/T^ = 0.1 when 
both the species have undergone Bose-Einstein conden- 
sation. Here is the condensation/critical tempera- 
ture under the mean-field approximation 'l3l for species 
1 with a total number of particles N . The mixtures pos- 
sess two distinct structures depending on the relationship 
between the external potentials on the two species. We 
have examined two cases, one with an equal energy spac- 
ing in the single-particle spectra and another with an 
identical external potential. We will elucidate the results 
for both cases. 

We first consider the case with the same energy spac- 
ing in the single-particle spectra, that is, wi = lo2- 
Two identical condensates are found when 7712/7711 ~ 1, 
N2/N1 — 1, and aii/a22 = 1, as shown in Fig. ^a). 
Comparing with the mean-field result on the one-species 
condensate , the double condensates in our simulation 
display the same expansion as Na.jk is increased from 2 
to 40.' 

In order to separate the two overlapped condensates, 
we perform simulations with various parameters 7772/7(7.1, 
Oil = 022, and N2 = Ni, respectively. Figures ^a)-(i) 
show the density profiles of the two species in the system 
of these cases. There are two significant features in the 
density profiles as 1712/1711 is increased from 1 to 4, as 
shown in Figs, ^a)-(c): The heavier species remains to 
be concentrated at the central area of the trap while the 
lighter particles are pushed outward and finally form an 
outer shell surrounding the heavier species. 

The interaction dependence of the density profiles is 
given in Figs, nd)-(f). The interaction strength can 
dramatically affect the particle distribution. Here we 
show the results for the same number of particles in 
each species with 7772/7771 = 4. The densities of both 
the species at the trap center decrease as the interaction 
strength is increased from = 0.15 to 0.3. However, 
under a stronger interaction strength, the lighter parti- 
cles are pushed out further away from the trap center. 
Heavier particles spread out with a less distance. 

Figures ^g)-(i) show the changes in the density dis- 
tribution as the total number of particles increases from 
200 to 1600 while keeping N2 = Ni, 1712/1711 = 4, and 
On = 022 = 0.15. The two condensates are always in 
partial overlap but with more the lighter particles going 



outward with a larger N. This means that when more 
particles are trapped in an experiment, the heavier par- 
ticles can be caged inside the light particles with both 
the species forming condensates at the low energy levels 
of the single-particle spectra. 

The condensation fraction in a species can be de- 
termined by counting the particles on the permutation 
rings ^3 • We have evaluated the condensation fractions 
for each species in the system with parameters shown in 
Fig.njand find them in the range of 0.8-0.95. This means 
that a reasonable fraction of particles is not in phase 
coherence with the rest in that species. Furthermore, 
particles can accumulate in any single-particle state, for 
example, the first excited state, to form a condensate. In 
fact, in certain parameter regions, the particles can be 
completely driven out of the single-particle ground state 
by strong correlation. 

To have a better picture of whether a particle is in 
its single-particle ground state or not, we further decom- 
pose the density of each species into that of the single- 
particle ground state and the corresponding excite states 
in Fig. 121 

We show here the percentages of the particles in each 
species in the single-particle ground states at different in- 
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FIG. 1: Density profiles of the double condensates with var- 
ious parameters at T/T^ — 0.1 and with A''i = N2. (a)-(c) 
show the changes of the density profiles at different mass ra- 
tios but with A'^ = 200 and an = 022 = 0.4. (d)-(e) show the 
cases at different interaction strengths with 7712/7711 — 4 and 
— 200. (g)-(i) show the dependence of the total particle 
number with m2/mi = 4 and an = 022 = 0.1. 
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FIG. 2: Decomposition of the density of each species into 
that of the single-particle ground state and the corresponding 
excited states at different interaction strengths with m2 /mi — 
4 and A''i = N2 — 100. The percentage in each plot indicates 
the fraction of the particles in the single-particle ground state. 



teraction strengths with 1112/1711 — 4 and Ni — N2 — 100. 
The percentage of the particles in the single-particle 
ground state drops as the interaction strength is in- 
creased for both the species. This can be interpreted 
as a result of competition between correlation and phase 
coherence. When the strength is small (with small hard- 
core radii), the particles are close to be free and they 
condense to the single-particle ground state when the 
temperature is lowered beyond the critical temperature. 
However, when the interaction strength is increased, par- 
ticles try to avoid each other to lower the total energy of 
the system, and more particles move into the low-lying 
excited states with higher single-particle energies. Our 
simulation shows that the lighter species is affected more 
drastically by the interaction strength than the heavier 
one because the corresponding percentage drops faster as 
the interaction strength is increased. 

Decomposition of the density profiles indicates that 
when the spatial separation between the two species oc- 
curs, the majority of particles from both the species do 
not occupy the single-particle ground states. This is be- 



cause when one species is pushed away from the central 
area of the trap by the correlation, the particles in both 
the species are also driven away from the single-particle 
ground states by the same force. The distribution of the 
particles in the energy space is strongly influenced by the 
energy spacing and the scattering lengths. Therefore, no 
significant separation between the two species can hap- 
pen in the energy space when ui = UJ2 and an = 022- 

So far we have shown that the lighter particles can 
form an outer shell surrounding the heavier particles. A 
mass ratio far away from 1, strong interaction, and large 
number of particles are all the contributing factors for 
the particles to accumulate in the higher single-particle 
energy states and to cause a spatial separation between 
the two species. However, the structure of a mixture 
turns out to be very different if the external trapping po- 
tentials on the two species are different from those with 
LOi = L02- Let us consider the case with an identical ex- 
ternal potential, that is, with k\ — 'm\iJ\ — k2 = m2w|. 
In Fig. 13 we compare the results under different exter- 
nal potentials. The percentage indicates the fraction of 
particles in the single-particle ground state. The simula- 
tion result shown in Fig. |2ta) under an identical external 
potential does not have any significant phase separation 
even though most particles are pushed out of the single- 
particle ground states. While keeping the same energy 
spacing, the spatial phase separation between the two 
species is obvious, as shown in Fig.|3Jb). 
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FIG. 3: Density profiles of two-species system with the iden- 
tical external potential and equal energy spacing. The per- 
centage in each plot indicates the fraction of particles in the 
single-particle ground state. A shell of lighter particles is 
formed under the equal energy spacing condition. 

More interestingly, the roles of the heavier and lighter 
species can be reversed under the identical potential case 
with the heavier particles forming an outer shell sur- 
rounding the lighter particles. To consider a possible 
experiment, we have performed the simulation with pa- 
rameters appropriate to the *^Rb-^^Na mixture. We set 
Na as species 1 and Rb as species 2. The total particle 
numbers iVi = N2 = 100 and the hard-core radii used in 
producing the results shown in Fig.^Ja) are equivalent to 
an actual mixture with s-wave scattering lengths of 3 nm 
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and 6 nm, respectively, and a total particle number on 
the order of 10^ 0. In Fig.^ the density profiles of both 
the species under the same single-particle energy spacing 
and identical external potential are shown. The energy 
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FIG. 4: Density profiles of Rb and Na witfi the identical ex- 
ternal potential and equal energy spacing. The percentage 
in each plot indicates the fraction of particles in the single- 
particle ground state. Separations in both the spatial and 
energy spaces occur with the identical external potential. 

spacing in the single-particle spectrum of Rb is about 
half of that of Na for the data shown in Fig. Efa) due to 
the difference in the trapping frequencies. The small en- 
ergy gap and large scattering length make it much easier 
for Rb atoms to be driven out of the trap center or its 
single-particle ground state whereas most Na atoms stay 
at the trap center or in its single-particle ground state. 
This is totally counterintuitive and in disagreement with 
the mean- field results No significant spatial sep- 



aration is found when the single-particle energy spacings 
are the same, as shown in Fig.^Jb); but increasing either 
the interaction strength or the total number of particles 
can result in the spatial separation with an outer shell of 
Na. 

Under the identical external potential, we have also 
found that the spatial separation of the two species is 
independent of the ratio TO2/TO1. Therefore, for a system 
containing any two species, it is energetically favorable to 
have the species with a larger scattering length to form 
a low-density outer shell surrounding the other with a 
smaller scattering length; this is consistent with the ex- 
perimental observation for the mixture of ^^Rb atoms 
with two distinct spin states 1] . It is also possible for a 
future experiment to realize the separation in the energy 
space as predicted here. 

Our simulations show that the external potentials can 
be used to tune the structure of a boson mixture. With 
the same trapping frequency for the two species, it is 
more likely for a system with a mass ratio close to 1 to 
form an outer shell of particles with a larger scattering 
length whereas for system with a mass ratio far away 
from 1 to form an outer shell of lighter particles. In 
the latter case, the large difference of external potentials 
between two species outweighs the moderate effect of the 
two-body interactions within a species. It is possible for 
the *^Rb-^^K mixture to form an outer shell of Rb atoms 
whereas for the ^^Rb-^'^Na mixture to form an outer shell 
of Na atoms and the future experiment will tell. 
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